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Abstract 

We consider the production of isolated prompt photons in hadronic collisions. We present a 
general discussion in QCD perturbation theory of the isolation criterion used by hadron collider 
experiments. The isolation criterion is implemented in a computer programme of the Monte 
Carlo type, which evaluates the production cross section at next-to-leading order accuracy in 
perturbative QCD. The calculation includes both the direct and the fragmentation components 
of the cross section, without any approximation of the dependence on the radius R of the 
isolation cone. We examine the scale dependence of the isolated cross section, the sensitivity 
of the cross section to the values of the isolation parameters, and we provide a quantitative 
comparison between the full R dependence and its small- R approximation. 
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1 Introduction 



The production of prompt photons with large transverse momenta p? at hadronic colliders has 
been the subject of a continuing effort, both experimentally JD, pi 0, S, M, pj and theoretically 
0, during the last fifteen years. The expression 'prompt photons' means that these 

photons do not come from the decay of hadrons, such as 7r°, 77, etc., produced at large transverse 
momenta. Prompt-photon production is an interesting observatory of the short-distance dynamics 
of quarks and gluons by using a hard colourless probe. It is complementary to the electroweak 
processes of Deep Inelastic Scattering and Drell-Yan pair production, and to pure QCD processes 
such as inclusive production of jets or heavy flavours. One of the main motivations of its study is 
its sensitivity to the gluon density inside the colliding hadrons [[u|. Indeed the gluon distribution 
function is involved already at the lowest order (LO) in the strong coupling a s through the 'QCD 
Compton' subprocess qg — > jq, which dominates at fixed-target energies. Moreover the point- 
like coupling of the photon to quarks in principle makes the process of prompt-photon production 
ideally free from the uncertainties inherent in jet reconstruction (as in the case of jet production) or 
in fragmentation of partons into hadrons (as in inclusive hadron production). 

In fact, the latter feature is not as ideal as initially thought, because prompt photons can be produced 
according to two possible mechanisms, one of them being a fragmentation mechanism. Whereas the 
contribution from fragmentation remains small (less than 10%) at fixed-target energies, it becomes 
dominant in inclusive prompt-photon production at colliders. This was already true at the CERN 
SppS, at least in the lower range of the px spectrum for the UA1 and UA2 experiments at the 
centre-of-mass (cm.) energy \fS — 630 GeV. It is true at the Tevatron (y/S ~ 2 TeV), and will 
remain true at the LHC (\/S = 14 TeV), since this dominance increases with \/S. 

To be precise, the collider experiments at the Tevatron and at the forthcoming LHC do not perform 
inclusive photon measurements. The background of secondary photons coming from the decays of 
7T°, 77, etc., overwhelms the signal by several orders of magnitude. To reject this background, the 
experimental selection of prompt photons requires isolation cuts. The isolation criterion used by 
collider experiments is schematically as follows. A photon is said to be isolated if, in a cone of radius 
R in rapidity and azimuthal angle around the photon direction, the amount of deposited hadronic 
transverse energy E^ ad is smaller than some value Ermax chosen by the experiment: 

E* ad <E Tmax inside (y - Vl f + {<j> - 7 ) 2 < R 2 . (1.1) 

In addition to the rejection of the background of secondary photons, the isolation cuts also affect 
the prompt-photon cross section itself, in particular by reducing the effect of fragmentation. 

Hadronic production of isolated prompt photons is not relevant only to QCD studies. The LHC 
experiments, ATLAS and CMS, will search neutral Higgs bosons in the mass range 80-140 GeV 
through the decay channel H — ► 77 |I3 . The production of pairs of isolated prompt photons with a 
large invariant mass is the so-called irreducible background to this search. Therefore, a quantitative 
understanding of this background is relevant to reliably estimate the expected significance S/y/~B 
(S and B being the number of expected Higgs boson events and background events, respectively). 
This provides another important motivation for studying isolation in prompt-photon reactions. 

Various issues on isolation of photons based on the above criterion]^] have already been discussed in 
the literature g |, 0, [T| ||. The aim of the present article is to provide a detailed QCD study 

2 Other isolation criteria have also been discussed. We mention, in particular, the alternative criterion proposed in 
ref. Jlq ] and studied in ref. Jig and in sect. 6.3 of ref. |L7[| . Th e related topic of isolated photons produced in e + e — 
collision s h as also been abundantly discussed. A variant of ( |L.l| ) suitable for e + e~ was studied in Jlj) and revisited in 
EU, Elf . A different criterion |22| for isolated photons in jets has also been applied in measurements of the LEP 
experiments. 
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of the isolation criterion (1.1), and to present a complete (without any small- R approximation) 
calculation at next-to-leading order (NLO) accuracy in QCD perturbation theory. 



The first theoretical question raised by the implementation of isolation is an issue of principle. Since 
the production of isolated photons, and more generally of any isolated particle, is no longer an 
inclusive process, one may question the validity of the factorization theorem of collinear singularities 
that is established in the inclusive case ]l9j . In the case of e + e~ annihilation into an isolated prompt 
photon plus hadrons, this issue was examined in ref. pi) . It was shown that factorization still holds 
in terms of the same fragmentation functions as appear in the inclusive case, whereas the dependence 
on the isolation parameters is consistently taken into account by the short-distance partonic cross 
section. The case of hadronic collisions is more involved, since factorization is complicated by the 
presence of initial-state collinear singularities in the perturbative calculation. In this paper we extend 



the treament of ref. pi) to hadronic collisions, we consider the isolation criterion (1.1) and some 



variants of it, and we discuss and prove factorization to any order in QCD perturbation theory. 
The evaluation of NLO QCD corrections to the inclusive production of a single isolated prompt pho- 



ton in hadron-hadron collisions was considered in refs. |9[ and 11 . The authors of ref. |S[] computed 
the NLO corrections to direct photon production, without including those to the fragmentation 
mechanism. The NLO calculation of ref. |ll[] includes both the direct and fragmentation mecha- 
nisms, and it uses the small-cone approximation R <C 1. More precisely, the approach of ref. |fTTf 
is based on an analytical treatment of isolation in which the transverse energies, rapidities and az- 



imuthal angles in eq. (1.1) are replaced by energies and relative angles defined in the hadronic cm. 
frame. This substitution of variables has no effect in the small-cone approximation. To go beyond 
the small-cone approximation, in this paper we carry on our study in terms of the variables 77, 4> 



and Et involved in the criterion (1.1). Since these variables are invariant under longitudinal boosts 
along the direction of the beam axis, the criterion in eq. ([O]) is current practice in experiments at 
high-energy hadron colliders. We perform a complete NLO calculation: it includes both the direct 
and fragmentation contributions without any small-cone approximation. Our approch is based on a 
combined analytical and numerical treatment, and it can be used to implement different variants of 
the isolation criterion ( |l.l| ). 

This paper is organized as follows. In sect. |^, we recall the two production mechanisms (direct and 
fragmentation) of a prompt photon. In sect. [| we present a brief pedagogical discussion on non- 
isolated prompt photons: using the collinear approximation, we derive analytic expressions for the 
direct and fragmentation contributions to the NLO cross section. The general factorization issues 
that are raised in the case of isolated particle hadroproduction are discussed in sect. 0. In sect. 0, 
we show how the implementation of isolation affects the analytic expressions of sect, pfin the limit 
of small values of the radius R of the isolation cone. To go beyond the small- R approximation, we 
implement the NLO calculation of the isolated-photon cross section in a computer programme of the 
Monte Carlo type. The computer programme, which enables us to evaluate the full R dependence 
of the cross section, is described in sect. [| where we also present some numerical results and 
comparisons with the small-i? approximation. In sect. |?j we summarize our results and we conclude 
with some remarks on questions left open when the cone radius R becomes small or the isolation 
condition becomes very tight (when ETmax becomes small). 



2 Mechanisms of prompt-photon production 

As mentioned in sect. |l|, schematically, a high-p-r prompt photon can be produced by two possible 
mechanisms: either it takes part directly in the hard subprocess, or it results from the collinear frag- 
mentation of a parton that is itself produced with a large transverse momentum. From a topological 
point of view, when a 'direct' photon is produced, it is most probable that it will be separated from 
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the hadronic environment, whereas a photon 'from fragmentation' is most probably accompanied by 
hadrons, except when the photon carries away most of the momentum of the fragmenting parton. 
These fragmentation configurations are rare and atypical, and they are precisely those that are not 
suppressed by the isolation criterion. 

The LO contribution to direct prompt-photon production is given by the Born-level processes qq — > 
jg and q (or q) g — ► "/q (or q). The computation of the NLO contributions yields 0(a s ) corrections 
coming from the subprocesses qq — > jgg, gq (or q) — ► jgq (or q) and from the virtual corrections to 
the Born-level processes. 

The calculation of these higher order corrections also yields the LO contribution of the fragmen- 
tation type (sometimes called 'bremsstrahlung contribution'), in which the photon comes from the 
collincar fragmentation of a hard parton produced in a short-distance subprocess (see, for instance, 
Diagram a). This contribution appears in the following way. A final-state quark-photon collinear 
singularity occurs in the calculation of the contribution from the subprocess gq — > 779. At higher 
orders, final-state multiple collinear singularities appear in any subprocess where a high-p^ parton 
(quark or gluon) undergoes a cascade of successive collinear splittings ending up with a quark-photon 
splitting. These singularities are factorized to all orders in a s and absorbed into quark and gluon 
fragmentation functions of the photon, D^(z, Mp) and D^(z, Mp), defined in a certain factorization 
scheme at a factorization scale Mp chosen to be of the order of the hard scale of the process. When 
the fragmentation scale Mp is large with respect to ~ 1 GeV, these functions behave roughly as 
a/a s (Mp), so that these contributions are of the same order in a s as the Born-level terms in the 
direct mechanism. Moreover, because of the high value of the gluon parton densities at small mo- 
mentum fraction x, the gq (or q) initiated contribution producing one photon from fragmentation 
even dominates the inclusive production rate in the lower range of the photon p^-spectrum at the 
Tevatron, and this will be even more true at the LHC. The calculation of the NLO corrections to 
the fragmentation contribution (see, for instance, Diagrams b and c) is thus required for a reliable 
and consistent treatment of prompt-photon production at this accuracy. 

In this paper, we call 'direct' the contribution given by the Born term plus the part of the higher- 
order corrections from which final-state collinear singularities have been subtracted according to the 
M S factorization scheme. The remaining part is called 'fragmentation' contribution, and involves 
the fragmentation function of a parton into a photon as defined in the MS factorization scheme. In 
our numerical calculations, we use the fragmentation functions of ref. p3[ . 
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Note, however, that the above distinction between the two mechanisms has no direct physical mean- 
ing beyond LO. From a theoretical point of view, the distinction is defined by an arbitrary choice. 
It follows from the necessity of factorizing the final-state collinear singularities and absorbing them 
into the fragmentation functions. This factorization requires the introduction of an arbitrary frag- 
mentation scale M F , which is an unphysical parameter. More generally, it relies on an arbitrary 
choice of the factorization scheme, which defines the finite part of the higher-order corrections that 
is absorbed in the fragmentation functions together with the singularities; the remaining finite part 
is then included in the higher-order contributions to the partonic cross sections. The dependence on 
this arbitrariness, and in particular on Mp, cancels only in the sum of the direct and fragmentation 
contributions, so only this sum is a physical observable. In particular, any experimental identifi- 
cation of direct and fragmentation contributions based on topological grounds does not match the 
'theoretical identification' used throughout this paper. 

3 Inclusive cross section 

We start the discussion by considering the inclusive cross section for the production of a non-isolated 
photon with momentum p 1 . The transverse momentum and the rapidity of the photon are denoted 
by Ptj and y 7 (since the photon is massless, y 7 = i] 7 where i] 7 is the pseudorapidity) , respectively. 
The hadronic cross section da / dpr -ydy 7 , which we denote by <j(p 7 ), is given by the sum of the 
'fragmentation' and 'direct' contributions. It can be written, in shorthand, as 

f 1 dz 

°(Pi)=Y, —v a (Py/z;n,M,M F )D2(z;M F )+a^(p 1 ; t i,M,M F ), (3.1) 

a J Z 

where a a and a 1 are the corresponding 'partonic' cross sections. The contribution a a describes the 
production of a parton a (a = q, q, g) in the hard collision, and DJ is the fragmentation function 
of the parton a into a photon. The direct contribution a 1 does not contain any fragmentation 
function; it corresponds to the point-like coupling of the large-pr photon to a quark produced in the 
hard subprocess. Note that a a and a 1 are not true partonic cross sections, since they include the 
convolution with the parton distributions of the colliding hadrons. 
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The cross sections ct 7 and a a are known up to NLO in a 

'a s (fi) 



^) vlornfo M) + (^) ^ho(PI M, M, M F ) , (3.2) 



d a (p;p,M,M F ) = 4™(p;M)+^j 4 (p; ft M,M f ) . (3.3) 

The expressions of a ho (and a Born) for the direct and fragmentation contributions can be found 
in rcfs. |t], |ll| and ref. ^4j, respectively. They depend on the renormalization scale /i, on the 
factorization scale M of the initial-state parton distributions and on the factorization scale Mp 
of the photon fragmentation function. These results have been the starting point of the prompt- 
photon phenomenology of the last fifteen years [Q . More recently, expressions involving resummation 



of logarithmic terms that are large at the phase-space boundary also became available |25|, 26 
§,[§[§■ 

To better understand the differences between the inclusive cross section and the isolated cross section, 
we need a more detailed expression of c(p 7 ) in eq. (3.1). Let us first consider the fragmentation 
contribution in the Born approximation: 



, brems, 1 1 f G i/A(x\, M) G j/B (x 2 ,M) 

STTb^ 'rrf J X\ Xo 



x D2(z,M F ) e( 

z z min) dz dr/i, drj a PTa dp Ta - (3.4) 

Here we are considering the photon-jet double inclusive cross section. The photon has transverse 
momentum pT-y = zpTa and rapidity r/ 7 = r/ a . The jet recoiling against the photon, here the parton 
b, has rapidity rjj et — 775 . |.M_b| 2 is the matrix element squared of the Born subprocess, averaged 
(summed) over initial-state (final-state) spins and colours. Gi/A{xi, M) and Gj/ B (x2, M) are the 
parton distribution functions of the incoming hadrons, with xi 2 given in the hadronic cm. frame 

by 

*i=^ (e^+e^) , (3.5) 

where transverse momenta, rapidities and the hadronic energy y/S are defined in the cm. frame. 
We rewrite (3.4) in a more compact form: 

x / E a — — dr) b drj a p Ta dp Ta , (3.6) 

J dp a dr/b 



where z m i n is given by 



Ell. ( e ±v-, + P ±v b 

Vs 



) . (3.7) 



We can write a similar expression for the higher-order (HO) corrections to the fragmentation con- 
tribution: 



da l S ms [A+B^ 1+j ^ 



X! J Ea d(JHO \^ +let+X ^ D'-[:.Mr) d: <///,,, <///„ //■; „ dp r „ ( .'■>.■ 



dp a drj 



jel 
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a,b,i Zmin z 



„ d(TBorn[A+B^a+b+x} 

x E a — — drjb drja p Ta dp T a (3.9) 

dp a drjb 



^ (^r) Y.\ dZ D 2(^ M F) K% ms [A + B^a+x](p Tl ^,M) . 



In eq. (3.E), the Mp dependence of the HO corrections has been made explicit. The remainder of 
the HO corrections is given by K^q 718 and no longer depends on Mp. The kernels Pi a (x) are the 
Altarelli-Parisi splitting functions. 

The direct contribution has the following form: 



da dlV [A+B->~f+jet+x] 



= 2ir 



+ 2tt 



dd B orn[A+B^+b] 



x E a 



dpj drj b 
dz e 2 a Pyq(z) In 

dd B or n 

\A+B->a+b\ 



drjb drj 7 pt 1 dpr- 



+ 2ir 



K HO 



dp a drjb 

[A+B^7+jet+x](p T •/, V, M) 



(PTj\ 

\Mp) 

drjb drja pr a dpr a 



(3.10) 



Here we have included QCD corrections up to NLO and also made the Mp dependence of the HO 
corrections explicit. The splitting function P J(I (z) of a quark q into a photon is 



P iq {z) 



l + (l-z) 2 



The LO evolution equations satisfied by the fragmentation functions are 



M 



dM F 



2vr 



2n 



( P «a + Pqa) ® D^(Mp) + P ga (g) D1{Mp) 



(3.11) 



(3-12) 



where the symbol ® denotes convolutions with respect to the momentum fraction z, and e a is the 
electric charge of the parton a. Using eq. ( 3.12| ), it is straightforward to verify the Mp -independence 
of the prompt-photon cross section up to NLO accuracy. 



A short description of the derivation of eq. ( 3.1C ) is instructive and helps to understand how the 
isolated cross section (see sect. ||) can be obtained from the inclusive cross section. Let us consider 
the subprocess l + 2^j + q + b, in which the photon is emitted by the final-state quark a (fig. |l|) . 
Its contribution to the cross section is 



da = 



dxi G\/ A (xi) dx 2 G 2 /b{x2) 



x ^\M\ 2 d Pl dp q dp b (2ir) n 6 {n) (pi +p 2 -p b -p-y -p q ) , 



(3.13) 



where we use the definition dp = d n p 5^ + \p 2 ) / (2i:) n 1 , and n = 4 — 2e is the number of space-time 
dimensions. We are mainly interested in the singular contribution (non-collinear contributions are 
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discussed in sect. ||) to cq. (3.13) coming from the configuration in which the final-state quark and 
the photon are collinear. Therefore the matrix element of the subprocess 1 + 2— > ~f + q + b can be 
approximated as 



|7W| 2 ~ 



(P yq (z) - ez) \Mi 

P 7 -Pq 

where Mb is the Born-level amplitude of the subprocess 1 + 2 - 
longitudinal variable z as follows: 

z = 



Pt- 



PTc 



(3.14) 

b in fig. [j], and we defined the 
(3.15) 



Inserting eq. (3.14) in eq. (3.13), we obtain 
V(n - 2)e 2 q fi 2e f dxx dx 2 



da 



2 3 (27r) 2 ™- 3 5' 2 



— — G{x x ,M) G(x2,M) (P 7 ,(z) - ez) \M B V Pr'-i d Pr 7 d V-y dq. 



( } \ d n - 2 Pq dj) q 



VS J \ J Pq ■ P~f 

where rji are the pseudorapidities, prb = \pt q +PTj\, and we have defined 



E 



PT-y 



e ±r >-< + PTq e ± ^ +PTbe ±Vb ■ 



(3.16) 



The factor V(n - 2) = 2tt : V /T (^ ) comes from the integration over the azimutha l ang le of the 
photon. The phase space in eq. ( p.lq ) can be written in terms of the variable z in eq. ( 3.15 ) and the 
difference <p q between the azimuthal angles of the quark and the photon: 



,2e 



Pq drj q 



P q -P~( 



V(n - 3)M 2e / ( sin( ^g) 2e d7 li p T q d PTg 

COSh(77 7 - T) q ) - COS (j) q pTqPT-j 



V>-3), 2 

\Pt~/ 



1 X ' r (sin(j> q ) 2e d(f> q dr/ q 



cosh(?7 7 — rj q ) — costftg \ z 



1 — z \ n 4 dz 



y 2 ■ 



(3.17) 



For the purpose of the discussion on isolated photons in sect. ^, we restrict the angular integration 
in eq. (3.17) to a cone of radius R in the ij-cj) space around the photon. This restriction does not 
change the singular collinear term (proportional to the pole 1 je) we are interested in. On the other 
hand, this produces a dependence on lni? 2 that we are going to use in sect. |5[ In the small-cone 
approximation defined by the constraint <d(R 2 — (f> 2 — (r/ 7 — f] q ) 2 ) with R 2 <§; 1, we can obtain a 
simple expression for the cross section. Using collinear expressions for the kinematic variables and 
defining the quark momentum p a — p-y + p q , we have 



PTa 
PTb 
>h 

Pi 



PT 1 

Z 

PTa , 

> 

PTj_ 

z 



+ e ±r ») 



(3.18) 
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Figure 1: Kinematics of HO corrections to the direct contribution. 




Figure 2: Kinematics of HO corrections to the fragmentation contribution. 
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The integrations over rj q and 4> q in ( |3.17 ) are easily performed, and we finally obtain (neglecting 
terms of O (R 2 )) 



da c 



V(n - 2)p n T -J fG 1/A (x u M) G 2/B (x 2 ,M) 
2 2 (2n) n ~ 2 S 2 



x {P iq {z) - ez) 



X\ X2 

i\ r(i-e) 



\M B \ P Ta d PTa dri a dr] b 
Anfj 2 ^ 6 



ej r(l-2e) \R 2 p 2 T (l-z) 2 



dz 



(3.19) 



where the variables Xi are given by 



x-l = Ell (e±»r + e**) 
2 zy/S V ; 



Equation (p . 1 Qj) can be rewritten as 



da c 



= 2tt 



a s (fj.) 



2 r l 



daBorn(A+B^a+b) 



dp a dr) b 



dr] b drj a pV Fa dp T a 



x {Dl(z,M F ) + —ei 



2|ii| E^-,) )p(:)| 



dz 



(3.20) 



where we have used the lowest-order definition of the quark into a photon fragmentation function 
in the MS factorization scheme: 



m[z,M F ,e) 



1 r(l-e) (\-K[i 



2\ e 



e T(l - 2s) V M| / 2tt 



-e q P^q(z) 



(3.21) 



We clearly recognize two parts in eq. ( 3.20 ) . The part proportional to D^{z, Mp), which comes from 
the collinear pole 1/e in eq. (3.19), is nothing but the Born-level fragmentation component in (3.4). 
The remaining part contributes to the direct component in eq. (3.1C). Since it has been calculated in 
the collinear approximation, it only contains terms that are either proportional to In R or constant 
in the limit R — > 0. Among the constant terms there is the term proportional to hx{pT-y/Mp), which 
was explicitly written in (3.1C). 



The contribution to the direct component of the inclusive cross section coming from the region 
outside the cone cannot be evaluated so straightforwardly, yet it must contain a term proportional 
to ln(l/i?) of the form 



da 



outside 



= 2vr 



Mm) 



2 r \ 



E a 



da B orn(A+B^a+b) 



dp a dr] b 



x dr] b dr] a p Fa d PTa-^e a 



dz 



(3.22) 



Indeed, the R dependence must cancel in the sum of the contributions from inside and from outside 
the cone. 

The NLO calculation described so far is extended to the case of isolated photons in sect. |^. Be- 
fore doing that, we discuss the validity of factorization of isolated photons at any order in QCD 
perturbation theory. 
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4 Isolated-particle production and factorization 



When no isolation criterion is applied, the inclusive photon cross section is computable by using the 
QCD factorization formula (3.1). The formula states that at large values of Ptj the short-distance 



dynamics is perturbatively computable in terms of partonic cross sections, while the dominant non- 
perturbative phenomena can be factorized in the parton densities of the colliding hadrons and the 
fragmentation function of the detected photon. Owing to the inclusiveness of the process, all the 
remaining non-perturbative effects have the form {Qq/pt^Y and are suppressed by some inverse 
power p of the photon transverse momentum. These terms are not indicated on the right-hand side 



of eq. (3.1), since they are negligible as long as PT-y is much larger than the typical hadronic scale 



Qo ~ 0(1 GeV). 



In the case of isolated photons, the isolation criterion (1.1) enforces additional phase-space restric- 
tions. This implies that the cross section is no longer fully inclusive over the hadronic final state 
and, hence, that the factorized expression fl3.l| ) is not necessarily valid |l(], [19|, |2(], Bl|. This issue 
was considered in ref . [ pl| , where factorization was proved at any order in QCD perturbation theory 
in the case of the isolation criterion of eq. ( |3.l| ). The discussion in ref. |2lJ was mainly focused 
on prompt-photon production in e + e~ collisions. To better clarify the differences between e + e~ 
collisions and hadron-hadron collisions, in the following we reconsider in more detail the latter case. 

In the theoretical literature (see e.g. refs. (h], 0) on isolated-photon production at hadron colliders, 



the isolation criterion (1.1) is sometime replaced by a similar one, in which the isolation is defined 
by an upper limit E max on the accompanying hadronic energy E had in the cm. frame rather than 
by an upper limit Exmax on the hadronic transverse energy E T ad . To discuss both criteria in a 
unified manner, we introduce a more exclusive definition of the isolation. The photon is isolated if 
the total hadronic four- momentum Q^f 1 inside the isolation cone is fixed at a given value Q^: 

Q h « d = Q» inside (y-y^f + <R 2 . (4.1) 

The cross section corresponding to this generalized isolation is denoted by da 9 ~ ls /d 4 Q. From it, we 



can recover the cross section o~ %s corresponding to the criterion in eq. (1.1) by simply integrating 
over as 

do®~ ^ s 

d A Q 0(E Tmax - Et) , (4.2) 

where Et denotes the transverse energy of the four- momentum Q^. We can also replace the hadronic 



transverse energy with the hadronic energy in the isolation criterion (IT), by simply replacing the 
transverse energies Et max, Et with the corr espo nding cm. energies E max ,E in the argument of 
the O-function on the right-hand side of eq. (|4.2| ). In the following we discuss factorization for the 
generalized isolated cross section da 9 ~ zs / 'd A Q. 

Our proof of factorization is based on the consistency, at the level of power unsuppressed con- 
tributions, between QCD perturbation theory and the full QCD theory. We exploit the formal 
correspondence between singularities in the calculation at the parton level and non-perturbative 
effects. Since the hadronic cross section is computable and finite in the full theory, the presence of 
singularities that are not cancelled or not factorizable at the parton level implies that of compensat- 
ing singularities of non-perturbative origin. Thus, we are led to consider the perturbative singular 
behaviour, which is due to the emission of soft and collinear partons. This picture obviously agrees 
with the physical expectation that factorization can only be spoiled by long-distance phenomena, 
such as those involving soft and collinear radiation. 

To give our formal proof, we consider the production of the prompt photon with momentum p 1 in 
the collisions of two hadrons Ha and Hb with momenta pa and ps- In general, the hadronic final 
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state contains n+l (with n > 0, since the transverse momentum of the photon has to be balanced by 
some final-state radiation) additional partons (hadrons) with momenta pi, ■ ■ ■ ,p n +i- Any isolation 
criterion applied to the photon is thus specified in terms of a function 



j-,(n+l) , \ 
F {ts} \PA,PB;PY,Pl,---Pn+l) 



(4.3) 



that depends on the momenta of the particles in the event. The subscript {is} denotes the depen- 
dence on the isolation parameters, whose precise definition is not explicitly spelled out for the mo- 
ment. The factorization of the prompt-photon cross section with isolation can be studied |3^, [}3| 
in terms of the properties of the isolation function in eq. (4.3). More precisely, factorization is valid 
provided F ^ s y fulfils the following requirements^] 



i) infrared safety: 

ii) collinear safety: 



(pa,pb;p 7 ;pi, ■ p i , - -Pn+i) o F$}(pa,pb;p*,;pi, ■ ■ -p n +i) , 



(4.4) 



F {is] {pa,pb;p~/-,pi, 



Pn+l) 



Pi II Pj 



Pi,pj, 

Hi) final-state collinear factorizability: 

F { ls} (pa,ps;p 7 ;pi, ■ ■ p t , ■ --Pn+i 

iv) initial-state collinear factorizability: 

F { is } {PA'Pb;PY>P1> ■ - Pi' ■ ■ -Pn+l 

Of course, the analogue of eq. (fh7|) with pa <-> 



F {is}(PA,PBiP-ylPlt 



,Pi +p j ,...p n+1 ) , (4.5) 



Pi II V~{ 



F {is} (PA,PB;p-r + pr,Pi, Pn+i ) , 



F 



(n) 
{is 



}(PA —Pi,PB',Py,Pl, Pn+l) ■ 



Pi II PA 

Pb is understood 



(4.6) 



(4.7) 



The requirement of i) infrared safety means that the cross section is insensitive to the momenta 
of arbitrarily soft particles. The requirement of ii) collinear safety implies that, when some final- 
state particles are produced collinearly, the cross section depends on their total m ome ntum rather 
than on the momentum of each of them. The general isolation criterion in eq. (4.1) depends on 
Qfi = Q'^ ad = J2j£coneP"j> tne total hadronic momentum inside the isolation cone. Since Q ^ d i s 
an infrared- and collinear-safe quantity, the isolation criterion in eq. (|4.l| ) fulfils eqs. (4.4) and (|4.5|) . 



The property in eq. (4.6) guarantees that all the long-distance phenomena related to the low- 
momentum fragmentation of the photon can be absorbed and factorized in the universal fragmen- 
tation function D2{z; Mp)- By universal we mean that it does not depend on the process and, in 
particular, it does not depend on the isolation parameters. Naive inspection of the isolation criterion 
in eq. (4.1) may suggest that it violates eq. (4.6). As pointed out in ref. plj , this is not the case. 

To explain the key point plj , let us first consider the following isolation criterion: 

n+l 



F 



(n+l) 



Ir}(Pa,Pb;PY,Pi, ■ --Pn+i) = <5 (4) Q c 



P~f 



Y,p 3 ®{R-Rj 



where 



R 



j- 



\J ivj - Vi) 2 + {4>j - <t>iY 



(4.8) 



(4.9) 



3 In eqs. (4.4)— bold-face characters are used just to emphasize the differences between left-hand and right-hand 
sides. 
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and R and Q^ Mt are external isolation parameters. The term in the square bracket on the right-hand 



side of eq. (4.8) is the total (hadron+photon) four-momentum inside the cone of radius R. Thus we 
have 

n+1 n+1 



Pi II p- 



(p 1 +p i ) + Y,Pj ®(R-Rj 



(4.10) 



3 = 1 



and the function F,^*) fulfils eq. (4.£) 



{Q"",-R} 



Then, we can observe that the momentum p 7 of the photon and the isolation parameter in 
eq. (4.1) are both kept fixed in the measurement of the cross section, so they can be regarded as 
external variables that are independent of the momenta PAiPBiP^iPx, ■ ■ -Pn+i- Therefore, by simply 
making the identification 



QT = Q, 



(4.11) 



the isolation criterion in eq. (4.1) can be recast in the form of eq. (4.S), which manifestly satisfies 
eq. (O). 



The properties in eqs. {iA)-{>Lt) are sufficient to prove factorization in the case of e + e~ collisions 
PH] , In hadron-hadron collisions, the cross section is affected by additional long-distance phenomena 
related to the non-perturbative binding of the colliding partons into the incoming hadrons. At the 
parton level, these phenomena lead to initial-state collinear singularities that have to be absorbed and 
factorized in the non-perturbative parton distributions of the hadrons Ha and Hb ■ The property 



in eq. (4.7) guarantees that the photon-isolation criterion does not spoil the factorization of the 
initial-state collinear singularities. Since the expression on the right-hand side of eq. (4.8) does not 
explicitly depend either on the incoming momenta pa,Pb or on any final-state momentum parallel 



to them, the property in eq. (4.7) is thus evidently fulfilled by the generalized isolation criterion in 
eq. (O). 



The main conclusion of our discussion on the generalized isolation criterion in eq. (4.1) is that QCD 
facto rization is valid at any order in perturbation theory. Factorization for the isolation criterion 
(1.1) with respect to the hadronic transverse energy (or to the hadronic energy) thus follows from 
eq. (fh^). Note, however, that da 9 ~ ls /d 4 Q fac toriz es as a fun ction of the fixed isolation parameter 
Q™' rather than as a function of (see eqs. (18) and ( 4.11 )). As discussed below, this functional 
dependence has influence upon the kinematical structure of the factorization formula for the isolated 
cross sections. 

We first discuss the isolation criterion (|Oj) with respect to the hadronic transverse energy. Since 
d(j9~ ls I 'd 4 Q factorizes at fixed Q™' = Q M + p AI7 , the constrained integration in eq. (4.2) leads to a 
dependence on the variable E™ 1 : 



= Et max + Et 7 — Et max + PT 7 



(4.12) 



The inclusive cross section is thus a function on the photon momentum p 7 and on the isolation 
parameters R and E™ 1 . It is convenient to define the variable 



PTj 

TP CUt 

LL/rp 



PT-y 



< 1 



■ PT' 



(4.13) 



The inclusive distribution da lB / dpT -ydy 7 with transverse-energy isolation, which we simply denote 
by cr is (p 7 ; z c ,R), fulfils a factorization formula analogous to eq. (3.1): 



a IS (p 7 ; z c , R)=J2 [ - a a ' is ^,R; (j,, M, M F ) D2(z; M F ) + d^ s ( Pj ; z c , R; M , 
j Q z \ z z / 



M, M F ) (4.14) 
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E 



; / - d a ' ls -,R; fJ>, M, M F ) D1(z; M f ) 0(z - z c ) 
J z \ z z / 



(4.15) 



The fragmentation function D2(z; Mp) is the same fragmentation function as appears in the non- 
isolated case. In particular, it does not depend on the isolation parameters. The dependence on the 
latter is fully embodied in the subprocess cross sections a a ' is and a 1,ts , which respectively give the 
fragmentation and direct contributions to the hadronic cross section. We recall that the subprocess 
cross sections a a ^ s and a 1,ts are obtained by convoluting the parton densities of the colliding hadrons 
with the cross sections d^- ls and cr^' IS of the partonic subprocesses i + j — > a + X and i+j^j + X. 
We have 



a a ' is (^;^,R;h,M,Mf) = V / d Xl dx 2 G i/A (xi, M) G j/B (x 2 , M) 



(4.16) 



and a similar formula relates ct 7,is to a 



Note that, according to eqs. (4.8) and (4.2), factorization holds at fixed i?™*. Since z c is obtained by 



rescaling iJ™ with the factor pr-y, the variable z c is a scaling variable with respect to factorization. 



In other words, the partonic cross section in eq. (4.15) depends onp 7 /z and z c /z. In particular, since 
z c is constrained to be z c < 1 from eqs. (4.S) and (4.2), this constraint propagates to a a,ls (p 1 / 'z; z c /z) 



as z c /z < 1. We made this condition explicit in eq. (4.15) 



In current experimental practice, Et, 
rameter Eh defined by 



is sometimes expressed in terms of the dimensionless pa- 

Et max 



PT- 



This parameter is related to our scaling variable z c by 

1 



1 + Sh 



(4.17) 



(4.18) 



The partonic cross section s in e q. Q4.15|) can be expanded as power series in a s analogously to the 
fully inclusive case in eqs. ( 4.19| ) and ( |3.3|) . Actually, at the Born level, it is straightforward to show 
fldj that a 1 ' 15 and a a - ls exactly coincide with the corresponding expression for the non-isolated case, 
apart from the overall constraint z c j z < 1 mentioned above. Up to NLO, we thus have 



a*> is (p;z c ,R;»,M,M F ) = 
a a ' ls (p; z c ,R;(jl,M,Mf) = 



a s (M) 



Borji 



(p; M) 



' Bor 



Q! S (M) 



a s (/x) 



a]i™(p;z c ,R;n,M,M F ) 



a Ho(P' z c ,R;n,M,M F ) 



(4.19) 
(4.20) 



In the following sections we compute the NLO terms a 



HO and 17 HO- 



We can now consider the variant of the isolation criterion (1.1) with respect to the hadronic cm. 
energy. We can straightforwardly follow the previous discussion on transverse-energy isolation. The 



transverse-energy isolation parameter E™ 1 in eq. ( 4.12 ) and the scaling variable z c in eq. ( 4. IE ) have 



to be respectively replaced by the energy isolation parameter E and the scaling variable 



En 



+ E y 



(4.21) 
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Denoting by a E %s the energy-isolatio n va riant of the inclusive distribution a ls , we still have a 
factorization formula analogous to eq. ( 4.15 ): it is sufficient to make the formal replacements 



cr 4S (p 7 ;z c ,i?) 

V z z 



R-is 



(py,z E ,R) 



^;^,R;ti,M,M F 
z z 



(4.23) 



In particular, the fragmentation function DJ(z; Mp) is not affected by the change of the isolation 
criterion. Note, however, that the strict formal correspondence between the two variants of the 
isolation criterion does not extend to the cross sections of the partonic subprocesses i + j ^ a + X 



and i + j — ► 7 + X. In fact, although the convolution structure of eq. ( |4.16| ) is still valid, we have 
to perform the following replacement: 



^a,is[ 



p Z c \ E ' ( P Z E E cut J7J cut 

z z J \ z z Xl VS x 2 VS 

(4.24) 

Unlike cr°' IS , its energy-isolation variant d1j E ~ ls does depend on E cut j x\\[S and E cut j x 2 \[S . This 
additional dependence follows from having defined isolation with respect to the hadronic energy 
in the cm. frame. Such a definition is not invariant under longitudinal boosts along the beam 
direction, thus leading to an entangled dependence on E cut and on the energies x\ VS/2,x 2 VS/2 
of the partonic beams. In other words, the dependence on the isolation parameter E cut is not 
kinematically factorized [ElJ from the dependence on the momentum fractions in the parton densities. 



Despite the formal correspondence in eq. (4.23), beyond the LO, the functional dependence of a %s 
on z E and R is not the same as the functional dependence of a ls on z c and R. However, the 
transverse-energy fraction J2j Erj/Ex-f and the cm. -energy fraction ^jEj/E^ coincide as long 
as all the partons (hadrons) j inside the isolation cone are either soft or collinear to the photon 
direction. This implies that the two variants of the isolation criterion can substantially differ only 
if the isolation cone contains at least one hard parton that is not collinear to the photon. Such a 
kinematical configuration is suppressed both in the (soft) limit z c — > 1 and in the (collinear) limit 
R — > 0. In either of these limits, the two variants of the isolation criterion perturbatively coincide: 
the order-by-order perturbative calculations of a ls and a E ~ ls differ by terms that are of 0(1 — z c ), 
when z c — > 1, and of 0(R 2 ), when R — > 0. 



The validity of factorization implies that the partonic cross sections in eqs. (4.19) and (4.2C ) are 
computable in QCD perturbation theory. Nonetheless, their fixed-order perturbative expansions are 
not always well behaved. This is the case, for instance, in the kinematical configurations of highly 
isolated photons (1 — z c 1) and of very small isolation cones (R 1). When 1 — z c <C 1, the fixed- 
order expansion contains large double-logarithmic contributions, (a s ln 2 (l — z c )) n , of soft origin. 
When J?< 1, the fixed-order expansion contains large single- logarithmic contributions, (a s \nR) n , 
of collinear origin. The effects of these logarithmic contributions at NLO are discussed in Sect. ^|. 

Another general source of misbehaviour in the fixed-order expansion of perturbatively computable 
observables is the possible presence of integrable logarithmic singularities at some 'critical' points 
away from the soft and collinear boundaries of the phase space |34|] . Such singularities occur |l9|, 
p0[ pl| in the spectrum of isolated prompt photons produced in e + e~ collisions, and are located 
at the critical point 2E~ f — z E \/S (in e + e~ collisions, the isolation is defined with respect to the 



cm. energies). The partonic cross section a^j ls (or afj E ls ) in eq. ( 4.24 ) has an analogous critical 
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point at 2Ex-y — z c \J x\XiS j coshy* (or 2Et^ — z® \J x\x?.S I coshy*), where y* — y 1 — In yj X1/X2 
is the photon rapidity in the partonic cm. frame. Unlike the cm. energy \/S in e + e~ collisions, 
the partonic cm. energy y/§ = y/xix^S is not fixed and depends on the momentum fractions x\ 
and X2 of the colliding partons. Any possible integrable singularities in the partonic cross sections 
a^ ls thus disappear after integration over x\ and xi (see eq. ( 4.16Q ). In conclusion, there are no 
critical points in the single-photon i nclus iv e cros s secti ons a ts (p 1 ; z c , R), d 1,ls {p; z c , R; [i, M, Mp) and 
d a > ts (p; z c , R; \i, M, Mp) (see eqs. ( 4.15 ), ( 4.19 ) and ( 4.20| )) in which we are interested throughout 
the present paper. Note, however, that critical points and related integrable singularities occur in 
less inclusive distributions of isolated prompt photons produced in hadron-hadron collisions (see e.g. 
ref. H). 

The photon isolation procedure actually i mplem ented by the DO Collaboration 0] does not exactly 
coincide with the isolation criterion in eq. (1.1). A cone of radius R around the photon direction is 
still considered, but the upper limit Exmax is enforced on the hadronic transverse energy inside an 
annular region of width AR (AR < R), rather than on the hadronic transverse energy inside the 
whole cone. More precisely, eq. fll.l| ) is replaced by the following: 



E T ad (R) - E T Ma (R - AR) < Et, 



jihad 



(4.25) 



where Ej , (i?) and Eft 1 (R — AR) are the hadronic transverse energies in the two cones of radius 
R and R — AR, respectively. This criterion is thus specified by the isolation function 



F {Et 1 J> a *,R,AR}(pA,PB-,p> y ;pi, ■ ■ -Pn+l) = 61 Et, 



n+l 

E 



E Tj Q(R- Rfo) 6(i? J7 -R+ AR) 



(4.26) 



It is straightforward to check that the properties in eqs. (4.4)-(4.7) are fulfilled by eq. (4.26), so 
the criterion (4.25) fulfils factorization. However, the kinematical structure of the corresponding 
factorized cross section, denoted by a D0 ~ ts (p 7 ; E Tmax , i?, AR), is different from that in eq. ( 4.15| ). 
We have: 



DO-is 



ETmax, R, AR) = y / <T a ' lS ( ] ETmax, R, AR; , 

„ Jo z v z 



p,M,Mp)Dl{z;Mp) 



(T 



7,D0-is 



(P 7 ; E T max,R, AR; fi, M, Mp) 



(4.27) 



The differences between eqs. ( gig ) and ( p~27l ) arise from the fact that the criterion ( |4.25| ) does 
not constrain the hadronic transverse energy collinear to the photon. Thus the isolation parameter 
ETmax is not rescaled by the photon momentum fraction z when going from the hadronic cross 



section o- D °- l;5 on the left-hand side of eq. (4.27) to the partonic cross section ^j a < D °-' LS on the right- 
hand side. Correspondingly, the isolation procedure does not set any absolute lower limit (such 
as z > z c in eq. ( 4.1 5| ) ) on the momentum fraction z in the photon fragmentation function. In 
particular, at the LO the isolated cross section (j 00 ^ 18 exactly coincides with the non-isolated cross 
section in eq. Q3.1| ). Of course, higher-order contributions to the isolated cross section are different 
from those to the non-isolated one, and tend to suppress the direct and fragmentation components. 

We conclude this section with two general observations. We have explicitly discussed isolated-photon 
production only in the case of hadron-hadron collisions. However, our discussion straightforwardly 
applies also to photon-hadron and photon-photon collisions: it is sufficient to substitute the parton 
densities of the colliding hadron for those of the colliding photon. Analogously, the final-state 
isolated photon can be replaced by any final-state isolated hadron (e.g. a pion) by substituting the 
fragmentation functions of the photon for those of the hadron. 
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5 Isolated cross section 



After the all-order proof of factorization, we explicitly study how the isolation criterion (1.1 ) modifies 



the inclusive cross section calculated at NLO. In particular, we show how the partonic cross sections 
&flQ and <t C hq depend on the cone size R. We also show how the constraint on the transverse energy 
inside the isolation cone can be translated into conditions on the integration range of the longitudinal 



variables in eqs. (3.9), (3.10) and (3.15). In our calculation we follow a procedure similar to that in 
rcfs. H 0: we start from the inclusive cross section and 
fulfil the isolation criterion. 



subtract' the contributions that do not 



5.1 Direct contribution with isolation 

Considering the 0(aa 2 s ) contribution to the direct component (see fig. [I]), when the momentum p q 
of the final-state quark is inside the cone around the photon, the isolation criterion (1.1) can be 
written as 



PTq < £h Pt 



7 ' 



or, equivalcntly, 



Pt 



7 



> 



1 



Zc. 



(5.1) 



PTq+PTj i- + S h 

Therefore, the contribution to the isolated cross section is obtained from that to the inclusive cross 



section by subtracting from the latter the part that violates the constraint (5.1) 



In terms of the result in eq. (3.20), obtained within the collinear approximation, eq. (5.1) is imple- 
mented by the con ditio n z > z c . Therefore the direct contribution to the isolated cross section is 
obtained from eq. (3.20) by simply subtracting the term 



i direct, sub 

da ljet = 2lT {—^- 



E 



E a - 



dxf 



Born 



dr]b dr] a pVj, a 3 dpr a 



dp a drj b 



Z7T 



2 In 



Rpr-y (l-z)\ l + (l-z) 



dz . (5.2) 



We recall that expression ( |5.2| ) has been obtained in the approximation R <C 1. All the terms of 
0(R 2n ) (n > 1) have been neglected and will be calculated as discussed in sect. |^. 

Note that the particular choice Mp = R pt 7 eliminates the independence from eq. (|]^) . The price 
to pay is the introduction of an additional independence in other terms of the cross section, notably 
through the ./lip-dependence of the fragmentation function. 

Note also that when the m ome ntum p q of the final-state quark in fig. [I] is outside the isolation 
cone, the isolation c riter ion (1.1) does not enforce any additional constraint. In particular, the HO 
contribution in eq. ( 3.22 ) is left u nchanged in the isolated case. As mentioned at the end of sect. ||, 
the In independence of eq. ( 3. 22] ) exact ly can cels the In independence of eq. (|3.20[ ). Since in the 
isolated case, we have to add eqs. ( |3.20| ) and ( 3.22 ) and then subtract eq. (|5.2]), the cancellation of 
the In i?-dependence does not occur anymore. When R <C 1, the NLO calculation of the isolated 
cross section is proportional to lnl/i?, and therefore it diverges to +00 when R — > 0. 

This divergence is unphysical. Its appearance at NLO simply means that, as soon as R is sufficiently 
small (see sect. O), the NLO calculation is not physically reliable. To improve the reliability of the 
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fixed-order perturbative expansion, higher-order contributions proportional to (a s lni?)™ have to be 
computed when R becomes very small. 

Moreover, the calculation in sect. || also tells us that non-perturbative contributions must be taken 



into account when R is very small. For example, the integral in eq. (3.17) is an integral over the 
virtuality p 2 of the parton a in fig. |l|. The integral of the perturbative matrix element is performed 
over the range < p 2 = 2p q .p 1 < R 2 p\^ (1 — z)/z. However, non-perturbative hadronization effects 
become dominant as soon as p 2 < Qq ~ 0(1 GeV 2 ). A strict perturbative treatment of the isolated 
cross section is thus justified only when 

R 2 P 2 Tl ^>Ql ■ (5.3) 

If R becomes too small, the perturbative calculation has to be supplemented by a careful treatment 
of non-perturbative phenomena. 

5.2 Fragmentation contribution with isolation 

To study the effects of isolation on the fragmentation component, let us consider the subprocess of 
fig. |^, in which a gluon with momentum p§ is emitted by the quark a that fragments into a photon 
with momentum p 7 . We denote by pj the momentum of the collinear hadronic fragments of the 
quark a. The gluon may or may not belong to the isolation cone around the photon, thus leading 
to two different restrictions on the allowed kinematics of gluon 5. 

i) The gluon is outside the cone. We only have a condition on p? f 

PTf<£hPT~t (5.4) 

Defining 

PTa = PT f + PTj , z=-^-, (5.5) 

PTa 



we write eq. (5.4) as 

z > z c . (5-6) 
The variable z is the fragm enta tion v aria ble appearing in the frag men tation function D^{z, Mf) 



Going back to expressions (|3.6|) and (|3.9| ), we see that condition (|5.4| ) restricts the z-integration 
range. The lower limit now is z c instead of z m i n (we are assuming that z c > z m i n , as is the case in 
the kinematical configurations of interest in experiments at high-energy colliders) . 



ii) The gluon belongs to the cone. The isolation criterion (1.1) implies 

p T f+PT5 < £hPr 7 , (5.7) 

or, cquivalcntly, 

PTa+PT5< — ■ (5.8) 

Zc 

In terms of the variable z in eq. ( |5.5| ) and of the variable x, 

PTa 



PTa +PT5 



(5.9) 



eq. (5.8) is written as 

z x > z c . (5.10) 
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The constraint in eq. (5.10) means that, when the gluon is inside the cone, we must also restrict the 
^-integration range, so that x > z c j ' z instead of x > z m i n /z has to be used in eq. ( |3.9| ). 

In summary, we obtain the fragmentation component with isolation by subtracting from the inclu- 



sive one the contribution that violates eqs. (5.6) and (5.10). Condition (5.6) is straightforwardly 
implemented by changing the z-integration range in eqs. (3.6) and (3.E). Condition (|5.10[ ) modifies 
the ^-integration range and the function K^Q ms in eq. (3.9). Working again within the small-i? 
approximation, we can perform a calculation similar to the one carried out in sect. [3] in the case 
of the direct component. We thus obtain the part of the higher-order correction that has to be 
subtracted from eq. (|3.9|) (after having implemented the constraint z > z c ) when R <C 1: 



do 



brems,sub 
j jet 



= 27T 



I D2(z,M F )dz 



J Z dxC F jln^ 



R PT-y (1 — x)\ / 1 



z M F 



1-x 



x E, 



do Bom (A+B->c+b+x) 

dp c dr) b 



drj b di] c p T c dp T < 



(5-11) 



Here the photon momentum is related to the momentum of the parton c by p 7 = zxp c . As in the 
case of the direct component, the explicit ^-dependence of the contribution in the curly bracket can 
be eliminated by choosing the scale Mp = RpTj, at the price of introducing independent effects 
through the Mi?-dependence of the fragmentation function DJ(z, Mp). 

In the above presentation, we have explicitly considered the subprocess in fig. ^, which is collinear- 
(when R — > 0) and infrared- (when z c — > 1) divergent and thus leads to the most important HO 
corrections. The same approach applies to other subprocesses. For in stanc e, parton a can be a gluon 
and parton 5 a quark. We obtain corrections similar to those in eq. ( |5.11 ), apart from the following 
replacement. The contribution inside the curly bracket on the right-hand side of eq. ( 5.11 ) is: 



In 



R PT-y (1 — X 



z Mp 



^)^)+~^r 4) (*) 



(5.12) 



where the terms 



p(n- 

qq 



A \x)=C F {l-x) 



build the unregularized Altarelli-Parisi kernel in n dimensions Pqq "* (x) 



Pqq 1 (x) 



HO corrections from the subprocess with generic parton species c and a in fig. 
eq. ( 5.11 ) through the substitution 



- £ P 9 ( r 4) (z). The 
are obtained from 



^ir 4) (z) 



The explicit expressions of Pat\x) and Pac can be found, for instance, in ref. J55|. 

The calculation described in this subsection is valid in the limit R <C 1 and thus neglects corrections 
of 0(R 2n ). These corrections are discussed in sect. ^|. 



(n-4), 



6 Numerical calculation and 0(R 2n ) terms 

The isolation criterion fll.l| ) modifies the calculation of the HO corrections to prompt-photon produc- 
tion by restricting the available phase space to final-state radiation. In sect, pi we have implemented 
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these cuts by working in the collinear approximation for pedagogical purposes. To go beyond the 
collincar approximation and keep the complete -R-dependence of the cross section, we have imple- 
mented all the LO and NLO contributions in a computer programme, according to a combined 
analytical and Monte Carlo approach. In this section we first briefly describe the programme, then 
we present the results of our numerical study. 



6.1 Brief presentation of the programme 

The code we use is derived from the NLO Monte Carlo programme DIPHOX |3f|, designed to 
calculate the double-inclusive cross sections 

da(A+B->Fi+F 2 +x) 
hi - 



dpi drj 2 

and associated distributions, where F\ and F2 are large-pr particles, photons or hadrons |35| p6|| . 
This programme combines the phase-space slicing method |3l|| and the subtraction method 
(3^, |33) to treat the soft and collinear singular parts of the perturbative matrix elements. 

For a generic parton subprocess 1 + 2— > 3 + 4 + 5, the photon and one outgoing parton, say 3 and 
4, have a high pr and are well separated in phase space, while the remaining final-state parton, 
say 5, can be either soft or collinear to one of the other four partons. The phase space is sliced by 
using two arbitrary, unphysical parameters pTm and Rrh, with prm <C \pt 3,4 1 and Rtii *C 1, in the 
following way: 

- Part I corresponds to \pr s| < PTm- This cylinder supports the soft and initial-state collinear 
singularities. It also yields a small fraction of the final-state collinear singularities from the 
subregion in which parton 5 is very soft. 

- Part II. a corresponds to the region where \pt b\ > PTm and the parton 5 is inside a cone C3 
about the direction of particle 3, defined by (y 5 — y 3 ) 2 + (</> 5 — 3 ) 2 < R\ Yl . This region 
supports the final-state collinear singularities appearing when 5 is collinear to 3. 

- Part II. fr is defined in a similar way as II. a, but with the replacement of particle 3 by particle 
4. The corresponding cone C4 supports the final-state collinear singularities appearing when 
5 is collinear to 4. 

- Part II. c is the remaining region: \px s| > PTm, and px 5 belongs to neither of the two cones 
C3, C4. This slice yields no divergences, and can thus be treated directly in four space-time 
dimensions. 

Collinear and soft singularities appear in parts I, II. a and II.fr. They are first regularized by dimen- 
sional continuation from 4 to n — 4 — 2e (e < 0) space-time dimensions. Then, the n-dimensional 
integration over the kinematic variables (transverse momentum, rapidity and azimuthal angles) of 
parton 5 is performed analytically over the phase-space regions in I, II. a and II.fr. After combina- 
tion with the corresponding virtual contributions, the soft singularities cancel, and the remaining 
collincar singularities that do not cancel are factorized and absorbed in the parton distributions and 
fragmentation functions. The resulting quantities correspond to pseudo cross sections, where the 
hard partons in the regions I, II. a and II.fr are unresolved from the soft or collinear parton 5, which 
has been 'integrated out' inclusively on these parts. The word 'pseudo' means that they are not 
genuine cross sections, as they are not positive-definite in general. These contributions as well as 
the one from region II. c are then encoded in the Monte Carlo computer programme. 

The integration over the phase-space region in II. c, which yields no divergences, is perfomed numer- 
ically without any approximation. The implementation of the isolation criterion is straightforward: 
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conditions (^^) or ( |5. 10 ) just cut the numerical phase-space integration when parton 5 is outside 



the isolation cone or inside the annulus R% h < (j/5 — j/ 7 ) + (^5 — 7 ) < i? 2 , respectively. 

The region in part I is treated according to the phase-space slicing method. The integration is carried 
out analytically by neglecting terms that are proportional to powers of prm and thus vanish when 
PTm 0. This approximation implies that, when using the numerical programme, the unphysical 
parameter prm has to be chosen sufficiently small for the results to be independent of pr m - In 
ref. | {35j] it was checked that pxm values of the order of half a per cent of the minimum of ptz and 
PT4 are sufficient. 

In principle, the dependence on the unphysical parameter R t h can be treated as that on PTm- 
However, such a procedure would lead to numerical instabilities when the radius R of the isolation 
cone is also small. Therefore, to avoid any approximation of the -Rt/j-dependence, the integration 
over the phase-space regions in part II. a and II. b is performed by using the subtraction method. The 
integrand (i.e. the square of the matrix element) is first replaced by its collinear approximation (see 
sect. ||), and the corresponding integration is performed analytically. This calculation is formally 
identical to the one presented in sects. |^ and |E|, the isolation-cone radius R now being replaced 
by the slicing parameter Rth- Therefore, after factorization of the collinear singularities in the 



fragmentation functions, the numerical programme contains expressions like (cf. eq. (5.11)) 

( a s (n) \ 3 sr f 1 rat i\ r \ ^ t? dd B orn{A+B^b+c+x) 

27r > / DJ{z,M F )dz E c — — di] b dr] c p T c dp T c 

\ n J t L„ dp c drj b 

« <* Br) (^) + + £7- + s" - ->} - ■ 

in which the isolation criterion is simply implemented by changing everywhere z m i n by z c . Then, 
to take into account the exact dependence on Rth, the difference between the exact matrix element 
and its collinear approximation is integrated numerically, yielding finite contributions (which are of 
0(R%%) when R t h — > 0) that are included in the Monte Carlo programme. The results, therefore, 
do not depend on R t h- 

In the following we present some numerical studies that implement the isolation criterion in 



cq. (1.1). However, our NLO Monte Carlo programme can also be used to implement other definitions 



of isolation, such as, for instance, the isolation criterion in eq. (4.25) 



6.2 Numerical results 

We start the presentation of our numerical investigations by studying the sensitivity of the cross 
section with respect to variations of the isolation parameters. We then examine the dependence 
of the NLO results on the factorization and renormalization scale. We also compare our complete 
numerical treatment with that in the collinear approximation. Finally, we study the isolated cross 
section as a function of the transverse momentum of the photon. 



6.2.1 Sensitivity to the isolation parameters 

We use exact NLO expressions for the isolated cross sections, where 'exact' means that the full 
independence is kept, i.e. all the terms proportional to R 2n {n > 1) are taken into account on top 
of the collinear approximation used in sect. ^. We choose kincmatical parameters corresponding to 
those used by the CDF experiment at Tevatron Run lb §: \/S = 1.8 TeV, -0.9 < y 1 < 0.9 and 
PT-f — 15 GeV, which corresponds to the lower range of the photon pr spectrum. We consider the 
prompt-photon inclusive cross section 
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da fVm-ax fa 

dy 



dpr 7 J Vmin dy 1 dp T , 



Similar studies can be done for photon-jet cross sections |38 . We use the NLO parton distribution 
functions of the set MRST-99 (3^] , and the NLO fragmentation functions of set II in Bourhis et al. 
]2l| . The calculations are done with Nf = 5 flavours. The renormalization and factorization scales 
fi and M are both set equal to j>t 7 /2. 

Table 1 shows the sensitivity of the cross section to the value R of the isolation cone. In this study 
we fixed £/,. = 2/15 ~ 0.13333, which means that all events with hadronic transverse energy larger 
than 2 GeV in the isolation cone are rejected. The results without isolation are also reported for 
comparison. We verify that the Born cross sections are not sensitive to the isolation radius, as they 
should. 



Isolation radius 


Direct contribution 


Fragmentation contribution 


Total 


R 


Born 


NLO 


Born 


NLO 


NLO 


1.0 


1764.6 


3318.4 


265.0 


446.7 


3765.1 


0.7 


1764.6 


3603.0 


265.0 


495.0 


4098.0 


0.4 


1764.6 


3968.9 


265.0 


555.6 


4524.5 


0.1 


1764.6 


4758.2 


265.0 


678.9 


5431.1 


Without isolation 


1764.6 


3341.1 


1724.3 


1876.8 


5217.9 



Table 1. Isolated cross sections (the values are given in pb/GeV) corresponding to Eh — 0.13333. 



It is interesting to note that the HO contributions, both to the direct and to the fragmentation 
components, increase when R decreases. This is due to the fact that the implementation of isola- 
tion amounts to subtracting a contribution proportional to In R from the non-isolated cross section 
(see eqs. (f[|) and (^TTTj) ) . Since this subtracted contribution is negative when R < 1, the HO 



contribution to the direct component of the isolated cross section is quite large for small values of 
R. A similar behaviour is observed in the HO contribution to the fragmentation component. When 
all contributions are taken into account, the total cross section (direct + fragmentation) strongly 
increases with decreasing R. 

In particular, when R = 0.1, the NLO calculation gives an unphysical result: the isolated cross 
section turns out to be larger than the non-isolated one! Such a behaviour had to be expected in 



view of the discussion at the end of sect. |5.l| . The NLO results in table 1 imply that the value R ~ 0.1 
is sufficiently small to demand the inclusion of beyond-NLO perturbative terms and non-perturbative 
contributions. 

The sensitivity of the cross sections to variations of Sh is displayed in table 2. Now we fix R — 0.7. 
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PiTlPTCTV flit 


Direct contribution 


Fragmentation contribution 


Total 


£h 


Born 


NLO 


Born 


NLO 


NLO 


0.03333 


1764.6 


3820.9 


60.3 


168.7 


3989.6 


0.06667 


1764.6 


3734.3 


135.2 


303.4 


4037.7 


0.13333 


1764.6 


3603.0 


265.0 


495.0 


4098.0 


0.33333 


1764.6 


3434.1 


571.2 


883.9 


4318.0 


0.66667 


1764.6 


3359.5 


930.4 


1307.5 


4667.0 


1.00000 


1764.6 


3340.4 


1173.1 


1579.4 


4919.8 



Table 2. Isolated cross sections (the values are given in pb/GeV) corresponding to R = 0.7. 



Note that, already at the Born level, the fragmentation component is quite sensitive to and 
strongly decreases when such a cut is installed. 

The ratios NLO/Born increase when Sh decreases, indicating that the effect of higher-order correc- 
tions is larger at small than at large £/,. This is due to the following mechanism. Radiation collincar 
to the photon is more suppressed by the transverse-energy isolation cut than hard non-collincar 
radiation. Since the collincar contributions are negative when evaluated in the MS factorization 



scheme (see e.g. eq. (3.20)), their strong suppression leads to a sizeable NLO correction when eu 
decreases. We also note that, as expected, the effect of the Eh cut-off is very large on the fragmenta- 
tion component, in which a large part of the z-integration domain is suppressed. The total (direct 
+ fragmentation) cross section at NLO is rather stable with respect to en variations, because of the 
behaviour of the direct contribution. 

We point out that there is no infrared divergence coming from the NLO collinear contribution 
when Eh — * (z c — > 1). The direct component contains a term proportional to jl dz ln(l — z) ~ 



(1 — z c ) ln(l — z c ) when z c — > 1 (see eqs. (3.2C) and ( p.2[ )). The fragmentation component involves 
terms proportional to (1 — z c ) ln 3 (l — z c ): a factor ln 2 (l — z c ) comes from integrating the behaviour 



ln(l —x)/(l — x) in eq. (5.11), and a factor (1 — z c )ln(l — z c ) comes from the convolution with the 
large- z behaviour, proportional to ln(l— z), of the NLO fragmentation function D^(z, Mp) pi}] . Still 
in both cases the NLO cross section is finite when Eh — > 0. However, there are infrared-divergent con- 
tributions that are not smoothed out by the convolution with the fragmentation function D^(z, Mp). 
They correspond to soft gluons emitted non-collinearly to the photon (for instance, emitted from 
the initial-state quark of the LO direct subprocess q + q — > 7 -I- <?) and produce terms proportional 
to R 2 In £h in the direct component of the NLO cross section |ll| . The effect of these terms is thus 
suppressed when the size of the isolation cone is relatively small. The numerical stability of the 
NLO results in table 3 suggests that the contribution of infrared-divergent terms is not dominant, 
unless the parameter Eh becomes very small. 

6.2.2 Scale dependence 

Until now all calculations have been performed with the renormalization scale fj,, the initial-state 
factorization scale M and final-state fragmentation scale Mp, all equal to PT7/2. Here we study the 
sensitivity of the isolated cross section with respect to scale variations. Choosing standard isolation 
parameters, Sh — 2/15 and R = 0.7, we vary the scale [i = M = Mp between PT7/2 and 2px-y- The 
results are given in table 3. They can be compared with the corresponding results obtained without 
isolation, which are displayed in table 4. 
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Scale 


Direct contribution 


Fragmentation contribution 


Total 




(NLO) 


(NLO) 


(NLO) 


p T /2 


3603.0 


495.0 


4098.0 


Pt 


3155.3 


576.2 


3731.5 


2p T 


2840.7 


631.7 


3472.4 


Table 3. Scale dependence of the isolated cross section in pb/GeV. 


Scale 


Direct contribution 


Fragmentation contribution 


Total 




(NLO) 


(NLO) 


(NLO) 


Pt/2 


3341.1 


1876.8 


5217.9 


Pt 


2643.7 


2122.6 


4766.3 


2p T 


2188.8 


2250.7 


4439.5 



Table 4. Scale dependence of the non-isolated cross section in pb/GeV. 



The cross section decreases monotonically in both cases when the scale increases. The magnitude 
of the effect is comparable in the two cases: a relative variation of the cross section by ~ 16% is 
observed in the range of scales considered. We conclude that the implementation of isolation does 
not induce any extra sensitivity to the scale arbitrariness, at least when the three scales fi, M, Mp 
are kept equal. 



6.2.3 Comparison with the collinear approximation 



The isolation criterion used in ref. ]TT[ is based on an upper limit on the hadronic cm. energy 
rather than on the hadronic transverse energy inside the isolation cone. As pointed out in ref. ]TT[ 
and discussed in general terms in sect. Q in the case of soft and collinear parton radiation, the two 
variants of the isolation criterion become equivalent by simply i dentif ying the cm. -energy parameter 
e of ref. [ fill with our transverse-energy parameter Eh in eq. ( 4.17 ). Since the approximate NLO 
calculation performed in ref. jy} is actually based on the soft and collinear limits, it can thus be 
compared with our calculation by setting Eh = £• 

More precisely, the authors of ref. |ll| derived approximate analytical expressions for the NLO terms 
cr^o in eqs. (4.19) and (4.2C). These expressions are valid in the small-cone^ approximation R <C 1, 
and thus neglect corrections that are of 0(R 2 ) when R — * 0. In this limit the collinear approximation 



is valid, and our results in eqs. ( |5.2| ) and ( 5.11[ ) fully agree with those in ref. p| . In the case of 
the direct contribution ctJ^q, the authors of ref. flllfl also computed the dominant correction due to 
soft-gluon emission. This correction behaves as Rrlneh when Eh — > and R<^1. 



R 


e h = 2/30 - 0.06667 


e h = 2/15 ~ 0.13333 




Small-cone approx. 


No approx. 


Small-cone approx. 


No approx. 


1.0 


292.5 


267.7 


482.4 


446.7 


.7 


313.7 


303.4 


512.8 


495.0 


.4 


344.0 


345.6 


560.9 


555.6 


.1 


431.8 


432.9 


678.9 


678.9 



Table 5. Comparison between our results (no approximation) and the small-cone approximation 
Jil"| for the fragmentation component. The cross sections are given in pb/GeV. 

4 In ref. [|n| the cone size <5 is defined with respect to the relative angles in the cm. frame. When ii< 1, we thus 
have R = 5 cosh r? 7 . 
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In the case of the fragmentation component, Table 5 shows a comparison between our NLO calcu- 
lation and that of ref. We see that there is a good agreement between the two calculations 
when R -C 1, the collinear approximation becoming more accurate as R decreases (our results are 
obtained by Monte Carlo integration with an accuracy of 1%). For large values of R ~ 1, the small- 
cone approximation overestimates the fragmentation contribution to the isolated cross section by 
about 8%. In the case of the direct component, similar results were found in ref. JTT[ | by a numerical 
comparison with the calculation performed by a Monte Carlo code. 

6.2.4 Effect of the isolation as a function of pt 7 

Finally we study the effect of isolation as a function of pr-y- Figure || shows the ratios 'isolated'/'non- 
isolated' of the direct contributions, the fragmentation contributions, and the total contributions to 
the NLO cross sections da/dpT 7 - The ratios are displayed in the range 15 GeV < px-y < 60 GeV, 
for two different choices of isolation parameters, both with an isolation cone of radius R = 0.7: 

one with a fixed value of E^ max — 2 GeV, 

one with a fixed value of Sh = Exmax/pT^ — 2/15 ~ 0.1333, 

which coincide when px 7 = 15 GeV. 

As for the fragmentation component, isolation using a fixed Exmax acts more and more severely 
as Ptj increases, while by fixing Eh the amount of accompanying transverse energy in the cone 
increases with ptj- Therefore the ratios 'isolated'/'non-isolated' with fixed Exmax and with fixed 
eh have opposite variations when px 7 varies. 

Isolation has a rather small effect on the direct contribution in both cases, since it does not act at 
the Born level. 

The effect of isolation on the total contribution to the NLO cross section slightly decreases as pr 7 
increases, but it actually depends only weakly onpr 7 . 

7 Conclusions and outlook 
7.1 Conclusions 

In this paper we have considered the production of isolated photons in hadronic collisions. We 
have shown that isolation does not spoil the factorization properties that are valid in the case of 
inclusive production, and we have discussed the general factorized structure of the isolated-photon 
cross section. We have then performed a detailed NLO study of the implementation of the isolation 
criterion in the direct and fragmentation components of the production cross section. For pedagogical 
purposes, we have first given analytic expressions of the isolated cross sections obtained by using the 
collinear approximation, i.e. valid in the limit of small radius R of the isolation cone. To go beyond 
the collinear approximation, we have implemented the isolation criterion in a computer programme 
that encodes the exact dependence of the cross section on the cone radius at NLO accuracy. 

The main result of our NLO numerical study is that the inclusive cross section for isolated photons 
has a magnitude comparable to the one for non-isolated photons. Each separate component, direct 
and fragmentation, sizeably depends on the isolation parameters, both at LO and NLO. However, 
the measurable cross section, given by the sum 'direct + fragmentation', does not vary by more 
than 10% when eh varies between 2/30 ~ 0.06667 and 0.333, or when R varies between 0.4 and 0.7. 
Nonetheless, with small cones such as R ~ 0.1, the result of the NLO calculation for the isolated cross 
section becomes larger than the one in the non-isolated case. This counterintuitive and unphysical 
result reflects the fact that the fixed-order perturbative calculation is no longer reliable when R is 
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Figure 3: Ratios 'isolated '/'non-isolated' for the direct contribution, the fragmentation contribution 
and the total contribution to the cross section da/dprj at NLO. 
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very small: when a s In 1/R 2 ~ 1, perturbative contributions beyond the NLO have to be taken into 
account; when RpT-y ~ 0(1 GeV), also non-perturbative contributions are demanded. The scale 
dependence of the cross sections is similar in the isolated and non-isolated cases: the cross sections 
decrease monotonically by ~ 16% when the common scale \i — M — Mp varies between j>t/2 and 
2pT- Comparing our NLO calculation with its small-cone approximation pi, we find that the latter 
is adequate at small R and overestimates the cross section by about 8 to 10% at large R values, 
R ~ 1. Finally, isolation is found to weakly affect the pT7"dependence of the cross section da/dpTy 
in the range 15 GeV < p Tl < 60 GeV. 



7.2 Open problems 

We conclude with some remarks on several questions left open by our NLO treatment when the cone 
radius R becomes too small, as well as when the parameter Et max is small and isolation becomes 
very tight. 



7.2.1 Small cone radius 

As seen numerically in sect. [| the result of the NLO calculation for the isolated cross section at 
PT-y — 15 GeV violates the physical constraint cr(with isol.) < er(no isol.) when R ~ 0.1. This value 
of R is smaller than the ones relevant to experimental practice. However, the value R ~ 0.1 is not 
much smaller than the value R = 0.4 that is used in the most recent measurements at the Tevatron 
H and that is going to be used at the LHC [[14| . Therefore, we should be concerned about the 
actual value of R below which the reliability of the fixed-order perturbative calculation breaks down. 

Note that pr-y = 15 GeV and R — 0.1 imply Rpr-y ~ 1 GeV and a s {p 2 ) In 1/R 2 ~ 1 (recall that 



fi = PT7/2 is used in the numerical results of table 1). Therefore, as discussed at the end of sect. 5.1 
the scale Rpr-y is close to the non-perturbative region and, at the same time, higher-order corrections 
proportional to (a s In 1/ R 2 ) n can be relevant. Summation of the logarithmic dependence on R to all 
perturbative orders, combined with a careful study of the border-line between perturbative and non- 
perturbative regions, has to be undertaken to improve our understanding of the small- R behaviour 
of isolated-photon prodution. Work in this direction is in progress, and the results will be reported 
elsewhere. 



7.2.2 Tight isolation 

Formally, the perturbative calculation of the cross section of isolated prompt photons is infrared- 
divergent in the limit en —> 0. Indeed, eh = would imply that the isolation cone about the photon 
would become an absolutely forbidden region of the phase space for gluon radiation, no matter how 
soft it is, thus spoiling the cancellation of infrared singularities between real and virtual soft-gluon 
contributions. As discussed in sect. [| at NLO this divergent behaviour shows up as a logarithmic 
term proportional to a s R 2 hiSh in the direct component of the cross section. 

In practice, using a cone of radius R — 0.7, we have found no significant infrared sensitivity in our 
numerical study, down to the very low value £h = 0.033333, which corresponds to Exmax = 0.5 GeV 



for a photon with pr-y — 15 GeV. Therefore the implementation of the isolation criterion (LI) with 
(finite but) tight transverse-energy cuts does not seem to destabilize the numerical convergence of 
the perturbative expansion. Nonetheless, owing to the presence of higher powers of lne^ at higher 
perturbative orders, the actual sensitivity of the cross section to very low values of Sh is probably 
underestimated in the present NLO calculation. 

In the present work we have treated all the quarks as being massless. This treatment is not adequate 
in the case of a heavy quark when Exmax is comparable to its mass. This issue especially concerns 
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heavy quarks that are experimentally not identified. In this case, better theoretical calculations and 
studies of heavy-quark fragmentation based on Monte Carlo event generators are required. 



In the NLO calculation, isolation is implemented at the parton level. In an actual event, a fraction 



of the hadronic transverse energy that accompanies a photon to be selected by the criterion (1.1) 
comes also from the low-p^ underlying event, as well as from pile- up effects in collisions at very 
high luminosity. Part of these effects can be mimicked in the NLO calculation by using an effective 
value of Exmax that is lower than that imposed at the detector level. However, this procedure 
does not take into account the fact that the underlying event causes also some suppression of the 
direct component at the Born level, which is independent of Exmax- Model estimates [|{J of this 
suppression indicate that it can be quantitatively relevant when the value of Et max at the detector 
level is small. In general, as the isolation becomes tight, the perturbative calculation has to be 
supplemented by careful studies of the effects of the underlying event and pile-up. 
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